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Abstract 

5h . 

We develop perturbative QCD factorization theorems for inclusive 6-hadron decays, in which radiative 
corrections characterized by the hadronic scale, the 6-hadron mass, and the W boson mass are absorbed into 
\Q " a heavy hadron distribution function, a hard b quark decay amplitude, and a "harder" function, respectively. 
Double logarithmic corrections associated with a light energetic final-state quark, which appear at kinematic 
end points, are absorbed into a jet function. Various large logarithms contained in the above functions are 
summed to all orders, leading to the evolution factors among the three characteristic scales. The heavy 
hadron distribution function is identical to the one constructed in the framework of heavy quark effective 
theory. It is shown that hadron kinematics must be employed in factorization theorems, and that perturbative 
contributions, depending on hadron kinematics, distinguish the lifetimes of the &-hadrons Bd, B s and Af,. 
Assuming the same heavy-quark-effective-theory parametesr Ai for these hadrons, we predict the lifetimes 
r(Bd) — 1.56 ps, t(B s ) = 1.46 ps and r(Ab) = 1.22 ps. We also predict the B u meson lifetime t(B u ) = 1.62 ps 
by varying the B meson distribution function slightly. All the above results are consistent with experimental 
^ I data. 
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I. INTRODUCTION 



<D 

' One of the puzzles in heavy hadron decays is the explanation of the low lifetime ratio r(Af,)/r(B £ i) = 
0.79 ± 0.06 0]. Inclusive heavy hadron decays involve both nonperturbative and perturbative corrections 
to the tree- level b quark decay amplitudes. In the framework of heavy quark effective theory (HQET) ||, 
nonperturbative corrections are expanded in powers of 1/Mf,, M& being the b quark mass, with the coefficient 
of each power proportional to a hadronic matrix element of local operators, and perturbative corrections are 
evaluated order by order at the quark level. The HQET prediction for the ratio r(Ab)/r(Sd) to 0(1/Mjj) 
is about 0.99 j|. When including the 0(l/Mjf) corrections, the ratio depends on six unknown parameters, 
and reduces to around 0.97 for various model estimates, which is still far beyond the experimental data. 
On the other hand, a phenomenological ansatz was proposed, in which the overall M| factor in front of 
nonleptonic decay widths is replaced by the corresponding hadron mass Mj^ Q]. This ansatz provides a 
solution to the puzzle, since (Mb/Ma,,) 5 — 0.73, Mb and Ma 6 being the B meson mass and the Af, baryon 
mass, respectively, is close to the observed ratio. As pointed out in 0], the same replacement also explains 
the absolute B meson decay rate, while the HQET prediction using the expansion parameter Mb accounts 
for only 80% of the decay rate. 

An alternative approach to inclusive heavy hadron decays is the perturbative QCD (PQCD) factorization 
theorem |(|. In this formalism radiative corrections to b quark decays are absorbed into different factors 
in factorization formulas for decay widths according to their characteristic scales. The soft corrections 
characterized by the hadronic scale of order Aqcd are factorized into a universal heavy hadron distribution 
function. The corrections characterized by the heavy hadron mass and by the W boson mass are factorized 
into a hard b quark decay amplitude and a "harder" function, respectively. The heavy hadron distribution 
function provides a summation over the initial states of the b quark, such that the virtual and real soft gluon 
corrections cancel exactly. The introduction of a distribution function indicates that a residual momentum 
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of the b quark, whcih arises as an effect from the light degrees of freedom in the heavy hadron, is allowed. 
The inclusion of the light degrees of freedom then demands the use of the heavy hadron kinematics, whose 
difference from the b quark kinematics sets a constraint on the magnitude of the residual momentum. 

The soft gluons are collected into the distribution function by employing an eikonal approximation, under 
which a b quark propagator is simplified into the propagator associated with the large component of the 
rescaled b quark field in HQET. This observation implies the identity of the heavy hadron distribution 
function to the one constructed in the HQET framework JtJ. Therefore, nonperturbative corrections also 
start from 0(1/ AI^) and the quark-hadron duality is respected in the PQCD formalism, the same as in HQET. 
However, as transforming the b quark kinematics to the heavy hadron kinematics through Mf, = Mh — A, the 
0{K/Mh) correction occurs M. It has been explicitly demonstrated, taking the B — > X U W as an example 
||, that the 0(1 /Mb) correction from the B meson distribution function cancels that from the phase space 
factor M%. 

As to perturbative corrections, various large logarithms contained in b quark decays are summed to 
all orders in the PQCD approach. The results are the renormalization-group (RG) evolutions among the 
three characteristic scales, such that the factorization formulas are independent of the renormalization scale 
/i. Double logarithms associated with light energetic final-state quarks are organized by the Collins-Soper 
resummation technique ||, leading to a Sudakov factor which smears end-point singularities and improves 
the applicability of PQCD. These perturbative factors, depending on kinematics of the initial heavy hadrons 
and of final states, differ among various decay modes of a hadron and among various hadrons. This is 
natural, because different decay modes should involve different energy releases. In the HQET approach 
perturbative corrections are evaluated to finite orders explicitly at the quark level with the renormalization 
scale fi set to a common value for all decay modes. It is then not a surprise that HQET predicts the lifetime 
ratio T(Ab)/r(-Bd) ~ 0.99, since nonperturbative corrections, supposed to distinguish the decays of different 
6-hadrons, are suppressed by l/M^ , and perturbative corrections to b quark decays are the same. 

The B meson distribution function has been extracted from the photon energy spectrum of the B — > X s "f 
decay ||, which provides the information of the b quark mass Mb and of the HQET parameter Ai for 
the Bd meson. Based on this information, we propose reasonable B s meson and Af, baryon distribution 
functions by assuming that they correspond to the same Ai. Convoluting these distribution functions with 
the perturbative factors, we predict the correct lifetimes of the 6-hadrons: r(Bd) = 1.56 ps, r(B s ) = 1.46 
ps and r(A h ) = 1.22 ps, and thus the correct ratios T(B s )/r(B d ) = 0.94 and T(A b )/r(B d ) = 0.78. By 
varying the B meson distribution function slightly, we obtain the B u meson lifetime t(B u ) = 1.62 ps, and 
the ratio T(B u )/r(Bd) — 1.04. It will be shown that perturbative contributions play an essential role in 
the explanation of the &-hadron lifetimes, and that our predictions are insensitive to the variation of the 
distribution functions, as expected from the quark-hadron duality. Our results of the semileptonic branching 
ratio £?sl = B(B — > Xlv) = 10.16% and of the average charm yield (n c ) = 1.17 per B decay are also 
consistent with experimental data. We present the semileptonic branching ratios and the charm yields in B s 
meson and A;, baryon decays as well, which can be tested in future experiments. 

In Sec. II the factorization theorem for semileptonic 6-hadron decays is constructed. The equivalence of 
the B meson distribution function to that derived in the HQET framework is demonstrated in Sec. III. 
The three-scale factorization theorem for nonleptonic decays is developed in Sec. IV. Various logarithmic 
corrections are summed to all orders using the resummation technique and RG equations in Sec. V. In Sec. VI 
we present the factorization formulas for the &-hadron decay widths and their numerical analysis. Section 
VII is the conclusion. The appendix contains the detailed derivation of the allowed phase space for heavy 
hadron decays. 

II. FACTORIZATION OF SEMILEPTONIC DECAYS 

We start with the factorization of the simplest case, the semileptonic decays B — > XJD, which correspond 
to the b — > civ decays at the quark level. These decays involve only the B meson distribution function 
and the hard b quark decay amplitudes due to the lack of the characteristic scale of the W boson mass. 
Nonperturbative dynamics is reflected by infrared poles in radiative corrections to quark-level amplitudes 
in perturbation theory. According to PQCD factorization theorems, these poles are absorbed into a hadron 
distribution function, which must be derived by nonperturbative methods or extracted from experimental 
data. A distribution function, being universal, is determined once for all, and then employed to make 
predictions for other processes containing the same hadron. In this section we demonstrate how to isolate 
infrared poles from radiative corrections and factorize them into the B meson distribution function. 
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Consider the one- loop corrections to the b — > civ decays shown in Fig. 1. Because both the b and c quarks 
are massive, there are no collinear (mass) divergences, and we concentrate only on soft divergences from 
vanishing loop momenta. The self-energy correction to the b quark in Fig. 1 (a) is written as 

v (a) ■ 2 n 2e f d 4 - 2 H 1 fa + M b fa- f+M b 2 2 

with Cf — 4/3 a color factor, p b and u b the b quark momentum and spinor, respectively, and p c = p b — q and 
M c the c quark momentum and mass, respectively, q being the lepton pair momentum. The S- function from 
the final-state cut on the outgoing c quark is included for the discussion of the soft cancellation between 
virtual and real corrections below. The b quark propagator (fa + M b )/(p b — M b ) after the self-energy 
correction helps the extraction of the soft pole in Eq. (l|). 

The soft divergence of E( a ) is isolated by the eikonal (I — > 0) approximation, 

vW -*n 2 e f d A - 2 H I fa + M b 2p» 

E softUfc = -tg Cf(i J j-^ T2 _ 2pt , [ + ^ _ M 2 Mfa, -M c ), (2) 

where the term / in the numerator and I 2 in the denominator of the b quark propagator have been neglected, 
and the equation of motion (jp b — M b )ub = has been applied to obtain the factor 2p b . We make the following 
expansion because of the on-shell condition p 2 — M 2 —* 0: 



1 1 Vl ~ M 2 

-2 Pb -l+p 2 -M 2 2 Pb -l 4( Pb -l) 2 



(3) 



where the first term on the right-hand side does not contribute, since it leads to an integrand with an odd 
power in I. The numerator of the second term removes the on-shell pole of the b quark propagator after the 
self-energy correction. Employing the equation of motion again, we obtain 



i 

(27r) 4 -~ e l 2 (v If 

where v = Pb/M b is the b quark velocity. Obviously, Eq. (|J) does not depend on the spin and mass of the b 
quark. 

The loop integral for Fig. 1(b) with a real gluon attaching the b quarks before and after the final-state cut 
is written as 

V (H i r 2e f d 4 ' 2 H A- /_ + Mb r fa- l + Mb fj, 

u b ^ U = -g C F , J J^-^ [pb _ l? _ M f {pb _ l? _ M f " 

x2nS(l 2 )S(( Pc - I) 2 - Ml) , (5) 

where T represents the other irrelevant vertices and propagators. The <5-function is associated with the final- 
state c quark with the momentum p c — I = p b — q — I in the case of real gluon emission. The soft divergence 
in Eq. (H) is also extracted by applying the eikonal approximation and the equation of motion to the two b 
quark propagators, leading to 

Z { 1 = -9 2 C F ^J ffl ^ 2^)5^ -If -Ml). (6) 

The above expression indicates that under the eikonal approximation, the b quark propagator is replaced 
by the eikonal propagator l/(v ■ I) and the quark-gluon vertex 7^ is replaced the eikonal vertex . These 
Feynman rules for an eikonal line are exactly the same as those associated with the large component h v (x) 
of the rescaled b quark field b v (x), 

1+ lb 

h v {x) = — - — b v (x) , b v (x) = exp(iM b v ■ x)b(x) , (7) 

that appears in HQET. This observation will become essential, as we demonstrate the identity of the B 
meson distribution functions constructed in the PQCD factorization theorems and in HQET. 
A straightforward calculation gives 
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41 = ^Cf-S(pI-M!), (8) 



The ultraviolet pole in S^ t has been subtracted, l/(— e) with e < denotes the soft pole, and the other 

finite terms irrelevant to our discussion have been dropped. The real gluon correction S^f t in fact contains a 
soft pole as p 2 c — > M^. Following the standard factorization of deep inelastic scattering, this pole is extracted 
through the convolution with a -B meson structure function. The introduction of such a structure function 
implies that the b quark momentum pb is allowed to vary. A variable pb is natural, since the light degrees of 
freedom in the B meson share various amount of meson momentum, such that the b quark is not always at 
rest. _ 

Define the residual momentum of the b quark as k — pb — MbV, v = (i> + , v" , vj_) = (1,1,0)/ V2 in the light- 
cone notation, for which k — (fc + ,0,O) is a convenient parametrization. The B meson structure function 
f(k + ) then describes the probability of finding a b quark with the residual momentum k + inside the B 
meson. Consequently, the one-loop corrections in Figs. 1(a) and 1(b) are modified into the convolutions of 
Eqs. ([l]) and (g) with f(k + ), respectively, where the c quark momentum is replaced by p c = MbV + k — q. It 
is legitimate to reexpress the factor l/(p 2 — Ml) 1+2e in Eq. (^|) as 

(pg - M*y+» = ^ 8{vl - Ml) + m -md + ■ (10) 

The second term in the above expression is defined via the convolution with an arbitrary function F(k + ), 

dk+ f ( fc+ ) _ ( d k+ F{k+) - m (11) 
in which 

- _ Ml - 2M b q° + q 2 - Ml 



2(M h v- - q- 



(12) 



is the value of k + determined by the on-shell conditionp^ = M%. Apparently, Eq. (|Tl|) is infrared finite as 
k + —>k {p 2 c — > Ml), and thus the second term in Eq. (10) is negligible. The first term in Eq. ([l(]) leads to 

V ( l = -^C F ±- e 5(pl-M?), (13) 

which cancels Sg"f t in Eq. (||). 

A similar cancellation occurs between the self-energy correction to the c quark in Fig. 1(c) and its corre- 
sponding real gluon correction in Fig. 1(d), when they are convoluted with the B meson structure function. 
The loop integrals are the same as those in Eqs. (|J) and (g) but with the b quark velocity v replaced by the 
c quark velocity v c = p c /M c . Using the eikonal approximation, the soft structure of the vertex correction in 
Fig. 1(e) with the virtual gluon attaching the b and c quarks is given by 

* - W / { t^ P :Z-i siA - M '> ■ (14) 

The soft structure extracted from the corresponding real correction in Fig. 1(f) is 

= 9 2 C F ^ J { C^v V lZl 2K6m(P ° ~ lf ~ K) ' (15) 

Following the similar reasoning, we show that and S^ft cancel exactly. The above soft cancellation 

can be generalized to all orders straightforwardly using the eikonal approximation, under which soft gluons 
detach from the quarks. In conclusion, the exact soft cancellation between virtual and real corrections to b 
quark decays must be implemented by introducing the B meson structure function. This is consistent with 
the Kinoshita-Lee-Nauenberg theorem, which states that infrared divergences in radiative corrections to a 
QCD process cancel when both final and initial states are summed. The B meson structure function f(k + ) 
simply provides a weighting factor for the summation over the initial b quark states. 
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The one-loop contributions to the hard b quark decay amplitude H for the semileptonic decays are defiined 
as the difference of the full radiative corrections and their eikonalized versions, that have been absorbed 
into the B meson structure function. As a demonstration, we perform the factorization of the self-energy 
correction to the b quark: 

(tree diagram) + Fig. 1(a) 
= (tree diagram) + Fig. 1(a) — (tree diagram) x E*^ + (tree diagram) x E^ t , 
= [(tree diagram) + Fig. 1(a) - (tree diagram) x E^jj t ] x (1 + E^ t ) + 0(a 2 s ) . (16) 

The first and second factors in the last line of the above expression belong to the first-order H and f(k + ), 
respectively. Extending the above procedures to all orders, we derive the factorization formula for the 
semileptonic decays, which is a convolution of the hard amplitude with the B meson structure function. 
Because of the soft subtraction (tree diagram) x E^., H is infrared finite and calculable at the quark level 
in perturbation theory. 



III. MOMENTS OF THE DISTRIBUTION FUNCTION 



A formal definition of the B meson structure function f(k + ) should reflect the soft dynamics associated 
with the B meson. Hence, we replace the b quark line by an eikonal line in the direction v, which collects 
infinite many soft gluons radiated by the b quark. From Eq. (Q) , the dependence on the b quark mass Mf, is 
removed by employing the rescaled b quark field b v . Specifying the large component h v of b v , the spin of the 
b quark remains fixed as radiating gluons, and its dependence also decouples. Therefore, the propagator and 
the gluon vertex associated with h v in HQET are exactly the same as the eikonal propagator and vertex, 
implying that f(k + ) is defined in terms of h v . The initial and final b quarks have a separation y~ in the 
minus coordinate, since the b quark momentum varies in the plus direction. At last, to render the definition 

gauge invariant, we insert a path-ordered exponential Pexp[i dtn ■ A(tn)\ in between the initial and final 
b quark fields, with n = (0, 1, 0) a vector on the light cone, forming 



/(fc + ) 



dy 
2tt 



(B(v)\h v (0)Pexp 



i I dtn ■ A(tn) 
>o 



h v {y-)\B{v)) 



(17) 



where \B(v)) is the initial state of the B meson. The Feynman rules for the path-ordered exponential are 
an eikonal line in the direction n with the propagator l/(n ■ I) and the vertex n M . 

The effects of the variation of the b quark momentum from the mass shell, determined by f(k + ), are 
nonperturbative. Because the variation is of order Aqcd> h is appropriate to expand the nonperturbative 
corrections contained in f(k + ) in terms of its moments [7j: 



dk+f(k+) = (B(v)\h v {0)h v (Q)\B(v)) = 1 , 
dk+k+f(k+) = (B(v)\h v (0)iD+h v (0)\B(v)) = , 
dk+k +2 f{k+) = (B v \h v (0)(iD+) 2 h v (0)\B(v)) = -Ax/6 



(18) 
(19) 
(20) 



where Eq . (|18| ) is the normalization of f(k + ), Eq. (|19| ) is the consequence of the equation of motion for h v , 
and Eq. (|20|) defines the HQET parameter Ai. Equation ( p^ ) implies that the 0(1/ Mb) nonperturbative 
correction to the semileptonic decays does not exist, as observed in HQET. That is, the quark-hadron duality, 
which states the equality of heavy hadron decay widths to heavy quark decay widths up to small difference 
of 0(1/ M 2 ) from Ai, is respected in the PQCD formalism. 

As argued before, the B meson structure function allows a variable b quark momentum, whose source is 
the effect of the light degrees of freedom in the B meson. This effect demands the consideration of B meson, 
instead of b quark, decays, and thus the B meson kinematics. Therefore, we transform the structure function 
f(k + ) into the usual B meson distribution function /s(z), z — /Pg, which describes the probability to 
find a b quark with the plus momentum zPg 0]. The relation between f(k + ) and /b(z) cab be extracted 
from their moments by employing the variable transformation, 
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Equation ( |l8| ) leads to the normalization of f B (z), 

f 1 M B 
j dzf B (z) = 1 , f B (z) = —j=f 



zM B - M b 



(22) 



Equation ( |l9| ) gives 



Mi - «)/bW = ( i - ^) /W) = ^ + ° (If J • ^ 



into which Eq. (^) has been inserted. Applying a similar manipulation to Eq. (|20|), we arrive at 

„ (1 . im „_^(i..^) +0 (^.). (24 , 

The HQET parameters A and Ai satisfy the relation 

A = M fl -M b + A-. (25) 

Therefore, the heavy hadron kinematics is introduced into the PQCD formalism rigorously via Eqs. (p"7|)- 
(pi|). It is then clear that the PQCD approach is equivalent to HQET in the treatment of nonperturbative 
corrections, and the moment A/ Mb is attributed to the replacement of the b quark mass by the B meson 
mass. An advantage of the B meson kinematics is that it provides the correct kinematic bounds. Take the 
lepton energy spectrum as an example. Using the b quark kinematics, the maximal lcpton energy Alb/2 is 
smaller than the correct value Mb/2. 

The B meson distribution function must be determined by means outside the PQCD regime. For its 
functional form, we propose 

z(l - z) 2 

f B {z) = N { j , (26) 

where the three parameters N, a and e are obtained from the best fit to experimental data. We have 
extracted Jb(z) from the photon energy spectrum of the radiative decay B — > X s "f, given by J^] 

0.02647 z(l - zf 

lB[Z) ~ \{z - 0.95) 2 + 0.0034z] 2 ' ( > 

The maximum of fs occurs at z ~ 1 as expected, since the b quark carries most of the B meson momentum. 
Substituting the above expression into Eqs. (22), ( p^ ) and (|24|), which relate the three parameters N, a and 
e to the HQET parameters, we obtain A = 0.65 GeV and Ai = -0.71 GeV 2 . 

As to the distribution functions for other 6-hadrons, we assume the parametrizations of the B s meson and 
A& baryon distribution functions, f Ba {z) and fi\ b {z), respectively, the same as Eq. (p6|). The moments of 
fB B (A b ) need to be treated as free parameters, since they have not been determined yet. The relations of the 
parameters N, a and e to the first three moments of fB B (A b ) ar e similar to Eqs. (p2|), ( |2^ ) and (pij): 

l 

fB s {A h ){z)dz = 1 , 

(1 - z)f Bsi A b) (z)dz = ^± +0 {Kl CTi /Ml s{H) ) , 

(1 - z) 2 f Bs (A b) (z)dz = 1 (a 2 Bs(A6) - iA 1Bs(Ab) j + 0(A 3 QCB /M Bs(Ab) ), (28) 



with 
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Aib,(a 6 ) - -6(B s (A i )\b(W ± yb\B a {A h )) 

M Bs{Kb) = M b + A Bs{Ab) - . (29) 

Note that only \ib b and AiA b are free parameters. Because of the parameters A = 0.65 GeV and Ai = 
-0.71 GeV 2 and M B = 5.279 GeV, the b quark mass M b = 4.551 GeV is fixed by Eq. (||). Substituting this 
M b into Eq. @ with M Bs = 5.369 GeV and M Ab = 5.621 GeV, A Bs and A A „ will be derived, if Ai Bs and 
Aia 6 are chosen. Using the values of A Bs (A b ) and X\ Bs f Ab ), combined with the normalization of f Ba (A b ), we 
obtain the corresponding parameters N , a and e. In the numerical analysis we shall assume that the HQET 
parameters Ai are the same for all 6-hadrons, i.e., X\ Bs = Ai Ab = Ai = —0.71 GeV 2 , which correspond to 
A B — 0.74 GeV and A Ab = 0.99 GeV. The B s meson and At, baryon distribution functions are then given 
by' 



0.02279z(l-z) 2 
JB ^ Z > [(z- 0.93) 2 + 0.0043z] 2 ' 



0.02095z(l- zf 

tKb{Z ' ~ [(z- 0.89) 2 + 0.0068z] 2 " ( } 

For the B u meson distribution function, the reasoning is a bit different. In fact, we can not distinguish 
the B u and Bd meson distribution functions, if f B is extracted from the photon energy spectrum of the 
B — » X s j decay, which acquires contributions from both mesons 0. To explain the B u meson lifetime in 
Sec. VI, we shall vary f B slightly: 

0.03358z(l-z) 2 

IBA ) [(z- 0.96) 2 + 0.0034z] 2 ' { 1 

which corresponds to A Bu = 0.65 GeV and Xi Bu = —0.81 GeV 2 , about 10% deviation from Ai for the B& 
meson. 



IV. FACTORIZATION OF NONLEPTONIC DECAYS 



The factorization theorem for nonleptonic B meson decays is more complicated than the semileptonic case. 
Nonleptonic decays involve three scales: the W boson mass Mjy, at which the matching condition of the 
effective Hamiltonian to the full Hamiltonian is defined, the characteristic scale t, which reflects the specific 
dynamics of a decay mode, and the b quark transverse momentum p±, which serves as a factorization scale. 
Below the factorization scale, dynamics is regarded as being purely nonperturbative, and absorbed into the 
B meson structure function, if it is soft, or into jet functions associated with light energetic final-state quarks, 
if it is collinear. Above the factorization scale, PQCD is reliable, and radiative corrections are absorbed into 
a hard b quark decay amplitude characterized by t, or into a harder function characterized by Mw- In this 
section we shall demonstrate how to construct the three-scale factorization theorem for nonleptonic B meson 
decays Q. 

Consider a nonleptonic b quark decay b — > cqq' through a W boson emission up to 0(a s ) virtual gluon 
corrections. The sum of the full diagrams does not possess ultraviolet divergences because of the current 
conservation and the presence of the W boson propagator. Reexpress a full diagram into two terms, where 
the first term is obtained by shrinking the W boson line into a point, and the second term is the difference 
of the full diagram and the first term. The first term corresponds to a local four-fermion operator (qq')(cb) 
appearing in the full or effective Hamiltonian with (qq 1 ) — 57^(1 — 75)9' the V — A current, and is absorbed 
into a decay amplitude H'(pb,Pj, /u), Pj being the outgoing quark momenta. H' is characterized by momenta 
smaller than the W boson mass Myy, since gluons in H' do not "see" the W boson. Note that the factorization 
in H' , which still contains the contributions characterized by the hadronic scale, is not complete yet. The 
second term, characterized by momenta of order M\v due to the subtraction term, is absorbed into the 
harder function H r (Mw, /*)• Following the steps that lead to Eq. (plf), we derive the factorization formula 
for the b quark decay width |llj , 

T = H r {M w ,fi) x H'fap^n) , (33) 

Though the full diagrams are ultraviolet finite, the factorization introduces the dependence on the renormal- 
ization scale n into H r and H', because the diagrams with the W boson line shrunk are divergent. 
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We then investigate infrared divergences from the diagrams contained in H'. These diagrams are basically 
similar to those in Fig. 1, where the four-fermion vertices have been adopted, and the final-state quark may 
represent the c, g, or q' quark. At this stage real gluon corrections are taken into account. If the final-state 
quark is c, the discussion the same as in Sec. II leads to the B meson structure function defined by Eq. Jl7|), 
which is thus universal. If the final-state quark is q, the loop integrals for Figs. 1(c) and 1(d) are written as 

1 1 \ n f d 4 ~ 2e l fo; — 7+771 h A + m 1 

(4 +m)E W = —ig Cpfj, J + m ) 7|1 

x5(p 2 -m 2 ), (34) 
+ / ^ 2 li i+ r 2 7,(A-- / + m)^ + "' 



x2tt5(; 2 )(5((^ -Z) 2 -m 2 ) , (35) 



with pj and m the final-state quark momentum and mass, respectively. 
Performing the loop integrations directly, we obtain 

a s _ / 1 3 , m 2 



E ' (C) = " | ^ =5 ' ^ " ™ 2) ' (36) 

(^ 2 P 2 ) £ 



* W =-^C F ,_^>>> . (37) 



Q 

7r ' (p 2 — to 2 ) 

The terms finite as to — > have been dropped and the ultraviolet pole has been subtracted in Eq. (|3 
Similarly, Eq. (J37|) is reexpressed, following Eq. (UQ), as 



E'W = -|^-U(p 2 -m 2 ). (38) 

Obviously, the soft poles e) disappear in the sum of Eqs. (^q) and ( |38| ) as expected. For m/0, there are 
no other infrared divergences, the same as in the c quark case. However, an additional infrared divergence, 
i.e., the collinear divergence, appears when the invariant mass p 2 of the final-state quark is equal to to 2 and 



vanishes, which remains even in the sum of Eqs. (36) and (p8[). This divergence comes from the region with 
the loop momentum / being parallel to pj and not small, for which both (pj — I) 2 and I 2 approach zero. Note 
that p 2 vanishes only in the kinematic end-point region. To absorb this end-point singularity, we introduce 
a jet function associated with the energetic outgoing light quark. Figures 1(c) and 1(d) are then trivially 
factorized into the jet function. 

For Figs. 1(e) and 1(f) with a gluon attaching the b quark and the q quark, the discussion is the same 
as for the c quark case, if the q quark is massive: there are no collinear divergences, and soft divergences 
cancel, as the two diagrams are convoluted with the B meson structure function. If the q quark is massless, 
the loop integrals for Figs. 1(e) and 1(f) are written as 

O = - V Cfc/, J — _r 4 _ l? _ M , V -Sip,) , (39) 

^,E'(/) = -g 2 C F ^J 'p-h^- mj^f^^6(l 2 )6((p 3 - I) 2 ) , (40) 

where T4 represents the four-fermion vertex. Similarly, collinear divergences from I parallel to pj with p 2 — > 

exist and remain after summing S and E (f\ These collinear divergences should be absorbed into the jet 
function associated with the q quark. 

The factorization of Figs. 1(e) and 1(f) in the collinear region requires an eikonal approximation. For 
convenience, we assume that pj is in the plus direction at the kinematic end point. When / is parallel to pj, 
only 7_ in the gamma matrices 7 M contributes, indicating that only 7~ in 7^ contributes, and that 7 in the 
numerator of the b quark propagator is negligible. Using the equation of motion for the b quark spinor, we 
have the eikonal approximation, 



Pb + Mj, „ 2p% ^ Pb 6 



n 



r = n ? , « - — 7 = ? • (41) 



(P6 - I) 2 — M 2 l 2 -2p b -l - Pb -l n-l 
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That is, the collinear gluon decouples from the b q uark, and attaches the eikonal line in the direction n 
denned in Sec. II. Using Eq. (|l|), Eqs. © and @ reduce to 



S™ = «W / ^SjpTfcCfe- /)^2^(( 2 )i(fc - . (43) 

Since the eikonal approximation in Eq. ( ptl] ) also holds for a massless quark propagator with M b — 0, the 
factorization of the other set of diagrams with a gluon attaching the q and q' quarks in the collinear region 
is similar: the collinear gluon decouples from the q' quark and attaches the eikonal line in the direction n. 
Extending the above analysis to all orders, the jet function is defined as the collection of infinite many gluon 
exchanges among the q quark and the eikonal lines in the direction n on both sides of the final-state cut. 
Note that the self-energy corrections to the eikonal line are excluded. The above conclusion applies to the 
other final-state q' quark: if the q' quark is massless, there exist collinear divergences at the kinematic end 
point, and another jet function needs to be introduced. 

With the jet functions and the B meson structure function defined in Eq. all the infrared divergences 
in nonleptonic B meson decays are extracted and factorized appropriately. The hard b quark decay amplitude 
H(t, fi) is then defined as the difference of H' and those diagrams which have been absorbed into the jet 
functions and the B meson structure funciton, and thus calculable in perturbation theory. Following Eq. (|l^) , 
we obtain the factorization of H', 

H'(p b , Pj ,fx) = H(t, (i) x J] Jjfan) x f(k+, M ) , (44) 

3 

where the index j runs over the light final-state quarks. Note that a trace of H and Jj is necessary, since 
the jet function carries the spin structure of the corresponding final-state quark. Combining Eqs. ( p3[ ) and 
(|44|), the three-scale factorization formula for nonleptonic B meson decay widths is written as 

T = H r (M w ,fi) x H(t,n) x Y[jj{pj,fJt) x f B {z,n) , (45) 

i 

where the structure function f{k + ,jj) has been transformed into the distribution function /b(z, /-*) 0- The 
fi dependence will disappear after performing a RG analysis, since a decay width is scale- independent. Note 
that the /x dependence of the B meson distribution function for the semileptonic decays was not considered. 
We shall show in the next section that such a /i dependence is not necessary in our formalism. 



V. LOGARITHMIC SUMMATIONS 



To regularize the collinear divergences, we associate extra transverse momentum p±, originating from the 
initial b quark, with the final-state quarks, which renders them off-shell by p±. This transverse momentum 
can be regarded as the factorization scale, above which RG evolutions are reliable, and below which dynamics 
is absorbed into the B meson distribution function and the jet functions. We work in the impact parameter 
b space, which is the Fourier conjugate variable of p±. Radiative corrections to nonleptonic decays then 
produce two types of large logarithms, hi(Mw/t) and ln(fb). Especially, the double logarithms hi 2 (pjb), pj 
being the large longitudinal component of the final-state quark momentum pj, appear in the jet functions at 
the kinematic end points. These logarithmic corrections should be summed to all orders using RG equations 
and the resummation technique J^] in order to improve perturbative expansions. The summation of the 
logarithms ln(Mw/£) is identified as the Wilson coefficients in the effective Hamiltonian, which describes the 
evolution from the characteristic scale M\y of the harder function to the characteristic scale t of the hard 
amplitude. The summation of the logarithms ln(tfr) leads to the evolution from t to the factorization scale 

The resummation of the double logarithms In (pjb) for a jet with a small invariant mass is written, in the 
b space, as f||l2f 

Jjip^b^fj,) = J J -(6,/*)exp[-2s(p_j,6)] . (46) 
The Sudakov exponent s is given by 
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Jl/b M 

where the factors A and B are expanded as 



ln(^)_Un.(//)) + UinJl,)) 



(47) 



AM=A M2l +A W(2lY , 



V J 3 7T V 2 



(48) 



in order to take into account the next-to-leading-logarithm summation. We adopt the one-loop running 
coupling constant, 



1 



ftln^/A 2 ) 



The above coefficients A^ and pi are 
33- 2n f 



Pi = 



12 



,m 4 4fa1 67 TT 2 10 8. . /e^ 

A (1 > = - , = n f + ~8 X In — 

3 ' 9 3 27 ; 3 f V 2 



(49) 



(50) 



where nt = 5 is the number of quark flavors, and 7 the Euler constant. The Sudakov factor exp(— 2s) 
exhibits strong suppression at large b, and approaches unity as pj < 1 jb. In this region the final-state quarks 
are regarded as being highly off-shell, and absorbed into the hard amplitude. Hence, double logarithms do 
not exist, i.e., exp(— 2s) — ► 1, and it is not necessary to introduce a jet function. For the detailed derivation 
of Eq. (46), refer to p2|Jl^ ]. The large scale pj will be chosen as the sum of the longitudinal components of 
Pj, i.e., pj = Pj + pj in the factorization formulas presented in the next section. 

The initial condition Jj{b, /i) of the resummation still contains single logarithms ln(6/i), which are summed 
by RG equations. The RG solution of the jet function is given by 



Jj (b, /x) = J exp 



-^7j(o<sO)) 
1/6 M 



(51) 



with the anomalous dimension 7^. It is more convenient to compute jj in the axial gauge n - A — 0, since the 
eikonal lines in the direction n disappear. The n dependence goes into the gluon propagator (—i/l 2 )N^(l), 
with 



n ■ I 



The lowest-order self-energy correction to the final-state quark is written as 



E 9 = -ig 2 C F [i 2e 



(2tt 



i4-2. 



-Ju. 



(pj-lf^P P 



(52) 



(53) 



which is obtained by replacing the gluon propagator —igV/l 2 in Eq. ([34]) by —iN^/l 2 . Extracting the 
ultraviolet poles from the above loop integral, we find that the term and the second term in Eq. ( |52] ) 
give — a s CV/(47r) x (1/e) and a s Cpl^ x (1/e), respectively. Their sum, S g = a s /ir x (1/e), leads to the 
quark anomalous dimension 7 g = —a s /ir in the axial gauge, jj is twice of 7 9 , i.e., jj = 2 r y q , because the 

self-energy corrections occur before and after the final-state cut. The initial condition jj , with the large 
logarithms collected by the Sudakov factor and by the RG evolution, can be approximated by its tree-level 
expression, i.e., the final-state cut. 

The RG solution of the B meson distribution function is 



/b(z>m) = /b 0*0 exp 



M ^ < <r\\ 
i/b A 1 



(54) 



where the initial condition /b( z ) absorbs the nonperturbative dynamics below the scale 1/b. The anomalous 
dimension 75 is also computed in the axial gauge, under which the path-ordered exponential in Eq. (|l7j) is 
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equal to unity. Hence, we need to consider only the self-energy correction to the h v field in Fig. 1(a). The 
loop integral is written as 



-ig C F n 



j4-2e 



I V„ 



(27r) 4 - 2e v ■ k v ■ (k - I) I 2 ' 
where the residual momentum k approaches zero. Using the expansion 



v ■ k 



v(k-l) v-l (v-l) 2 ' 



Eq. (p5|) reduces to 



ig CfH 



(2n)^J^l) 2 ~l 2 ~ 



(55) 



(56) 



(57) 



which is the same as in Eq. (^) but with the gluon propagator —ig^jl 2 replaced by —iN^ /I 2 . There- 
fore, Eq. (17) indeed generates the factor 1 + in Eq. (jlj) at first order. It is easy to find that the 
contribution from the second term in Eq. fl5^ ) vanishes, and that the g^ v term gives £„ = o. s Cf / {2w) x (1/e). 
Adding the self-energy corrections on both sides of the final-state cut, we derive the anomalous dimension 
7s = —a s CF/Tr of the B meson distribution function. 
The RG solution of H r is 



H r (M w , fi) = H r (M w , M w ) exp 



Mw 



-p-jHri&sifJ')) 



(58) 



with jh t the anomalous dimension of H r . The initial condition H r (Myy , Mw) can be safely approximated 
by its lowest-order expression Hr' = 1, since the large logarithms ln(Mvp//x) have been organized into the 
exponential, which is identified as the Wilson coefficient c(/z), 



c(p) = exp 



-irjHriotaQj.)) 



(59) 



For the explicit expressions of the Wilson coefficients, refer to JTJ] . 

At last, the anomalous dimension of the hard amplitude H is given by jh — —lH r — Y]j 7j — 7s because 
of the scale invariance of a decay width in the full theory. Applying the RG analysis, we obtain 



H{t,n) = H{t,t)exp< - 



* dp 



7H,.(a s (/i)) + ^2jj(a s (p)) +7s(a s (/x)) 



(60) 



The scale t will be chosen as the maximal relevant scales, 



t = max [ pj ; , — 



(61) 



Since the large b region is Sudakov suppressed by exp(— 2s) as stated before, that is, t remains as a hard scale, 
and the large logarithms ln(t//i) in H have been grouped into the exponential, the initial condition H(t,t) 
is calculable in perturbation theory. Hence, the Sudakov factor, though important only in the end-point 
region, improves the applicability of PQCD to inclusive nonleptonic heavy hadron decays. 

Substituting Eqs. @, (§l|), (§|), @ and @ into Eq. @, we derive the RG improved factorization 
formula for the nonleptonic B meson decay widths, 



T = c(t)H(t,t)f B (z)cxp\ -J2 



2s(p j ,b) 



dp 
i/b M 



7j 



dp, 
i/b fi 



7s (a* (m)) 



(62) 



where the cancellation of the \x dependences among the the convolution factors is explicit, and the two-stage 
evolutions from 1/6 to t and from t to Mw have been established. Note that Eq. (^) in fact denotes the 
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integrand appearing in the factorization formula, and that the variables t and b will be integrated out. The 
initial conditions jj ^ in Eq. have been incorporated into the evaluation of H(t, t). We emphasize that 
the Wilson coefficient appears as a convolution factor in the three-scale factorization theorem, instead of 
a constant (once its argument p is set to a common value for all decay modes) in the HQET approach. 
Since the perturbative factors depend on the hadron kinematics, they vary in different decay modes of a 
hadron and in different hadrons. ft will be shown that these perturbative effects play an essential role in the 
explanation of the 6-hadron lifetimes. 

In the semileptonic case there are not the Wilson coefficient c(t) and the Sudakov factor exp(— 2s) because 
of the absence of the double logarithms. Therefore, the large scales pj do not exist, and t is equal to 1/6, 
for which the RG evolution factors governed by the anomalous dimensions "fj and 75 disappear. This is the 
reason we did not perform the RG analysis for the semileptonic decays in Sec. II. 



VI. FACTORIZATION FORMULAS AND NUMERICAL ANALYSIS 



After developing the factorization theorems, we present the factorization formulas for the B meson semilep- 
tonic and nonleptonic decay widths. Define the momentum of the light degrees of freedom as p — (p + , 0, p_i_). 
The B meson is at rest with the mometum Pg = Mb/ V2(l, 1,0). The b quark momentum is then written 
as pb = Pb — P = (zAIb/V%, Mb/V%, — Pi), where the momentum fraction z = p^ '/Pg = 1 — V%P + /Mb is 
the same as that defined in Sec. III. The lepton and neutrino momenta involved in the semileptonic decays 
B(Pb) —> X c + I (pi) + v(p v ) are expressed, in terms of light-cone coordinates, as 



Pi = (pf,Pi ,0_l) , p u = (pt,P v ,Pui) , 



(63) 



where the minus component p l vanishes for a massless lepton. For convenience, we adopt the scaling 
variables, 



2E[ 
M B 



Q 



yo 



2<Zo 
M B 



(64) 



with the kinematic ranges, 



x ■ 



2yfa < x< 1 + a — (3 , 

ct < y<a+(l + a — (3 — x)- 

2(y -a) 



x + va? 



4a 



x — v? 



4a 



< Vo< 1 + V - P 



x + v 'x 2 — 4a 

where Ei is the lepton energy and q = pi + p v the lepton pair momentum. The constants a and (3 are 



(65) 



a 



Mf 
Ml 



ft 



__Ml 
-Ml ' 



(66) 



Mi and Mb being the lepton mass and the D meson mass, respectively. Mb arises as the minimal invariant 
mass of the decay product X c . For the derivation of Eq. (f35|), refer to the Appendix. 

The factorization formula for the semileptonic decay width is given, in the b space, by |15 



£sl 

r 



Ml 

2tt 



dxdydyo 



dz 



dbbf B (z)J c (x, y, yo, z, b)H(x, y, y , z) , 



(67) 



with To = (G f /16tt 3 )\ V c b\ 2 M B . The momentum fraction z approaches 1 as the b quark carries the whole B 
meson momentum in the plus direction. The minimum of z, determined by the condition p 2 . > M 2 , is 



V 

2 y 



Ml 
Ml 



V x 2 — 4a 



.y° + u. + - 

2 X X 



l _ m _ 
2 



\/x 2 



4a 



_yo y « 

2 x x. 



(68) 
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The function J c denotes the Fourier transformation of the final-state cut on the c quark line, which is in fact 
part of the hard amplitude. J c and the lowet-order H in momentum space arc fig] 



J c = 8{pl - Ml) , 



H = {pb -Pu)(pi -Pc) , 



x(l — z) 
V 1 x 1 — 4a 




Z X X 



(1-z) 



M 2 C - pi 



(69) 



(70) 



The universal B meson distribution function fg, determined from the B — * JG7 decay, has been given in 
Eq. @, which minimizes the model dependence of our predictions. 

For the nonlcptonic decays, we consider the modes b — > ccs and b — > cud. Ignoring the penguin operators, 
the effective Hamiltonian for the b — > ccs decay is 



H, 



efit 



4gj 

7f 



^K* s [ci(/i)Oi(M) + c 2 (/i)C2(/x); 



(71) 



with Gf the Fermi coupling constant, V's the Cabibbo-Kabayashi-Maskawa (CKM) matrix elements, the 
four-quark operators 0\ — (sb)(cc) and O2 = (c6)(sc), and the initial conditions c\(Mw) = 1 and C2(Mw) — 
0. For the b — * cud decay, V cs and the c and s quark fields are replaced by V u d and the u and d quark 
fileds, respectively. It is simpler to work with the operators 0± = ^(02 ± Oi) and their corresponding 
coefficients c±(/i) = C2(fi) ± ci(/x), since they are multiplicatively renormalized. In the leading logarithmic 





approximation c± are given by 



c±(m) = 



a s (M w ) 



(72) 



with the constants 2j + = —7- = —2 and = 5. 

To simplify the analysis, we route the transverse momentum p± of the b quark through the outgoing c 
quark as in the semileptonic case, such that the c quark momentum is the same as before. For kinematics, 
we make the correspondence with the c (u) quark carrying the momentum of the massive (light) lepton r 
(e and /z) and with the s and d quarks carrying the momentum of v. The scaling variables are then defined 
exactly by Eq. (Q) . The factorization formula for the nonleptonic decay widths is written, according to the 
three-scale factorization theorem in Sec. IV, as ll5f| 



1 NL 

"r7 



Ml 
2tt 



dxdydyo 



dz 



bdb 



N c + 1 2 



<(t) 



N r — 1 



c 2 -(t) 



xfB(z)Jc(x, y, 1/0, z, b)H(x, y, y , z)S(pj,t, b) 



(73) 



The B meson distribution function being universal, is the same as that for the semileptonic decays. The 
factor S is the result of the Sudakov resummation and the RG evolutions, given by 



S(pj,t,b) = exp< 



t ^ 

2s(pj,b) + I —jj(a s (fl)) 

l/b M 



1/6 M 



(74) 



with j = s for the 6 — > ccs mode and j = u, d for the b — > cud mode. The explicit expressions of pj are 



Ps 



Pd 



Mb 
V2 



(yo - x) , p v 



xM B 

~7T 



(75) 



The anomalous dimensions 7j = — 2a s /7T and 75 = —Cfc< s /it have been computed in Sec. V. It has been 
found that the single-logarithm evolutions governed by 7, and 75 enhance the nonleptonic branching ratios 
and thus lower the semileptonic branching ratios fl5fl , consistent with the observation in [[Hit . 
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The above formalism can be generalized to the B s meson and A;, baryon decays straightforwardly, for 
which the functional forms of the hard amplitudes H, the final-state cut J c , the Wilson coefficients c±, and 
the Sudakov evolution factor S remain the same, since they are evaluated at the quark level. The only 
differences arise from the replacement of the Bd meson mass Mb by the B s meson mass Mb, and by the A;, 
baryon mass M\ b , and from the heavy hadron distribution functions given in Eqs. (P) and @. We then 
proceed with the numerical analysis of the factorization formulas in Eqs. ( |67| ) and (|73|) for the semileptonic 
and nonleptonic 6-hadron decays, respectively, choosing the CKM matrix elements \V CS \ — \V u d\ = 1.0 and 
\Vcb\ = 0.044, and the masses M c = 1.6 GeV, M D = 1.869 GeV and M T = 1.7771 GeV. Our predictions are 
not sensitive to the change of M c : the difference of the results for M c — 1.5 GeV from those for M c = 1.6 
GeV is less than 5%. We obtain the lifetimes t{Bj) — 1.56 ps, t(B s ) = 1.46 ps and t(A(,) = 1.22 ps, or the 
ratios T(B s )/r(B d ) = 0.94 and r(A 6 )/r(B rf ) = 0.78. 

The experimental data of the 6-hadron lifetimes are summarized below. The lifetimes of the Bd and B s 
mesons and of the A fc baryon are r(B d ) = (1.55±0.03) ps, t{B s ) = (1.47±0.06) ps and r(A 6 ) = (1.23±0.05) 
ps, respectively, from the CERN e + e~ collider LEP measurements [Q. Recent CDF results yield t(A&) = 
(1.32 ±0.16) ps @. The lifetime ratios are then T(B s )/r{B d ) = 0.95±0.06 and r(A 6 )/r(fl d ) = 0.79 ± 0.06 
from LEP, and r(A&) /r(Bd) — 0.85 ±0.11 from CDF. Obviously, our predictions are almost the same as 
the central values of the LEP data. Employing the B u meson distribution function in Eq. ( |32| ) , we obtain 
t(B u ) = 1-62 ps or the ratio T(B u )/T(B d ) = 1.04, which is in agreement with the experimental data 
t(B u ) = 1.65 ± 0.04 ps in (k]]. 

To test the sensitivity of our predictions to the variation of \ib 3 and Aia,,, we adopt the HQET parameters 

A mcson _ ^baryon = _ Q 4 ^2 from q CD gum for ^ ^ ^ correS p 0ndiong \b 3 = 0.77 

GeV and A Ab = 1.03 GeV). The lifetimes t(B s ) = 1.49 ps and t(A 6 ) = 1.26, or the ratios T(B s )/T(B d ) = 0.96 
and T(Ab)/T(Bd) = 0.81, are derived, which are still much smaller than those from the HQET approach. 
This test indicates that our results are insensitive to the variation of the distribution functions (less than 
4% under 40% variation of the HQET parameters), and that the explanation of the experimental data is 
mainly due to the PQCD effects. To confirm this observation, we set all the 6-hadron masses in the phase 
space factors To and in the perturbative factors to the b quark mass Mb in Eqs. (|67| ) and (|73]) . That is, the 
difference among the factorization formulas for 6-hadron decays resides only in the distribution functions. 
The lifetime ratios T(B s )/r(B d ) — 1.01 and r(A(,)/r(B,j) = 1.09, which are even greater than unity, are 
obtained. Therefore, the perturbative contributions are indeed responsible for the low &-hadron lifetime 
ratios. 

The above results are aummarized in Table I. In Table I we also present the branching ratio of each mode in 
the Bd and B s meson decays and in the Af, baryon decays in terms of the quantities r TV = B(b — > ctv) / B{b — > 
cW), r u d = B{b — > cud)/B(b — > cW), r cs = BR(b ccs)/B(b — > dv), the semileptonic branching ratio i?sL 
and the average charm yield (n c ) per 6-hadron decay. It is observed that our predictions of B$l = 10.16% 
and (n c ) = 1.17 for the B meson are consistent with the experimental data: _Z?sl = (10.19 ± 0.37) and 
(n c ) = (1.12 ± 0.05) from the CLEO group @), and B S l = (11.12 ± 0.20) and (n c ) = (1.20 ± 0.07) from 
the LEP measurements It is also observed that B$l, (n c ) and all r's except r ud increase a bit with the 
6-hadron masses. It is interesting to examine these tendencies in future experiments. 

VII. CONCLUSION 

In this paper we have developed the PQCD factorization theorems for inclusive 6-hadron decays by care- 
fully analyzing the infrared divergences in radiative corrections to b quark decays. Radiative corrections 
characterized by the hadronic scale are absorbed into the heavy hadron distribution function or into the jet 
functions. Above the hadronic scale, perturbative contributions characterized by the 6-hadron mass and the 
W boson mass are absorbed into the hard b quark decay amplitude and the harder function, respectively. 
Various large logarithmic corrections have been summed to all orders using the resummation technique 
and the RG equations, leading to the evolution factors among the three characteristic scales, such that the 
factorization formulas are /^-independent. We have shown that the soft cancellation between virtual and 
real corrections demands the introduction of the heavy hadron distribution function, and thus the 6-hadron 
kinematics into our formalism. Using the PQCD factorization theorems, we have been able to explain the 
lifetimes of the 6-hadrons Bd, B u , B s and Af,. It has been confirmed that our predictions are insensitive 
to the variation of the distribution functions as expected from the quark-hadron duality, and that the low 
6-hadron lifetime ratios are mainly attributed to the perturbative contributions. 
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APPENDIX 

In this appendix we derive the phase space constraints for the decay b — > civ. The constraints for the 
decay b — ► ccs are the same with the c and s quarks corresponding to the I and leptons, respectively. The 
constraints for the decay b — > cud is the massless lepton case of the results presented below. 

We derive the constraint of x = 2Ei/Mh from the energy and momentum conservations associated with 
the decay i? — > X c /j>, 

M H =E x +Ei+E v , 

= px + P; + Vv ■ (76) 

Because of Ef = |p/| 2 + Mf , the minimum of £7/ occurs at the minimum of |p/|, |p/| = 0, and thus the 
minimum of x is given by 



_ 2Mj_ 
Xmin ~ M H ■ 

To obtain the maximum of x, we find the maximum of |p;| under Eq. ([76] 



(77) 



M H = V(P/ + P„) 2 + Ml + ^| P( | 2 + Mf + |p„| 2 , (78) 

into which px — — (p T + P^) has been inserted. It is easy to observe that the maximum of |p;| corresponds 
to p„ = 0, leading to 

(M 2 H + Mf - M 2 X ) 2 - AM 2 H Mf 

IPimaxI = 7772 ' ( 79 J 



and 



4M£ 



M 2 M| 



= 1 + T^-XFT- (80) 



Combining Eqs. @ and ©, we obtain 



2Mi Mf M 2 

irr- <x<l + ic7t-T7 



<x<l + ^--£, (81) 



where Mx has been replaced by the D meson mass Md, the mass of the lightest charmed hadron. 
We then derive the constraint of the kinematic variable y = q 2 /Mfj with 

q 2 = (Pi+Pv) 2 = Mf + 2(E t - \ Pl \co80)E u , (82) 

where 8 is the angle between the vectors p; and p„. The minimum of q 2 occurs, as cos 6* = 1 and E v takes 
its minimal value under the constraint from Eq. (|7q), 



M H -Ei = ^\pi\ 2 + El + 2\ Pl \E v cos6 + M 2 X + E u . (83) 

Obviously, the minimum of E u is zero, which can be achieved by increasing Mx- We then have 

(q 2 Un = Mf . (84) 

On the other hand, the maximum of q 2 occurs, as cos 9 = — 1 and E v takes its maximal value under Eq. (|8^), 
which corresponds to the minimum of Mx- Setting Mx = Md, we obtain 

Mfj + Mf -M 2 D - 2M H E l 

W " 2(M H - El - ^Ef-Mf) ' (85) 
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and 

(<Z 2 ) max = Mf + {Ml + Mf - M 2 D - 2M H E l ) El + ^ Et ~J Il \ . (86) 

M H - E t - y/Ef - Mf 

Combining Eqs. (|^) and (|S6]), we have 

M l , .. / M l d , M l M l 



Ml~ y ~ Ml{ 1+ Ml Ml X )\2- X -^lMfjMl 'J' ^ 

At last, we derive the range of the variable yo = 2qo/Mn with 

y = E l +E u , (88) 

or equivalently, the range of E v under the constraints of Eqs. ( |S2"|) and (jS3"l). Since E v decreases with cos 6 
in order that q 2 maintains constant, the minimum of E v occurs at cos# = — 1. Equation ( |82] ) then gives 

Eumin = mwEf^W)' m 

and 

Vo -I 2 fa- M ' 2 / M ^ (90) 

It is easy to observet that increases as Mx decreases. Setting Mx = Mr>, and substituting the expression 
of cos# from Eq. ( |8^ ) into (|83]), we obtain 

_ MI-2MhE 1 -MI+^ 

f-Jy max n * x ' \ / 



2M 



and 



M 2 

y max = l + y- T #. (92) 



Combining Eqs. (|90|) and (|92|), we derive 



2(y-M 2 /M 2 H ) _ ^ , M 



i + v/a; 2 - 4MJ/MI ~ M Z H 



<yv<i + y-T^ ■ (93) 
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Table I. The 6-hadron lifetimes r{Bd) 1 t(B u ), t(B s ), and r(Ah), and their ratios to T{Bd). Predictions of 
r Tv, r ud, r cs, B SL , and (n c ) for the 6-hadron decays are also listed. 







r TV r ud r cs B SL (n c ) 


T (ps) 


T/T(B d ) 


Ai = 


-0.71 GcV^ 


0.224 5.98 1.64 10.16 1.17 


1.56 


1.0 


AlB„ 


= -0.81 GeV^ 


0.222 6.11 1.54 10.14 1.16 


1.62 


1.04 


AlB s 


= -0.4 GeV^ 


0.231 5.25 1.44 11.21 1.16 


1.49 


0.96 


AlB s 


= -0.71 GcV^ 


0.237 5.30 1.69 10.85 1.18 


1.46 


0.94 


AlA b 


= -0.4 GcV^ 


0.254 5.14 1.75 10.94 1.19 


1.26 


0.81 


AlA b 


= -0.71 GcV^ 


0.261 5.26 1.68 10.87 1.18 


1.22 


0.78 
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